Mon. Not. R. Astron. Soc. 000, [HS] (2006) Printed 26 January 2013 (MN I^TeX style file v2.2) 



Gravity or turbulence? II. Evolving column density PDFs 
in molecular clouds 

Javier Ballesteros-Paredes^*, Enrique Vazquez-Semadeni^, Adriana Gazol^, Lee W. 
Hartmann^, Fabian Heitsch^, and Pedro Colin^ 

^ Centra de Radioastronomia y Astrofisica, Universidad Nacional Autonoma de Mexico, 
Apdo. Postal 72-3 (Xangari), Morelia, Michocdn 58089, Mexico 



2 



Department of Astronomy, University of Michigan, 500 Church Street, Ann Arbor, MI 48105, USA 



^ Department of Physics and Astronomy, University of North Carolina Chapel Hill, 
CB 3255, Phillips Hall, Chapel Hill, NC 27599, USA 



Submitted to MNRAS, 26 January 2013 



ABSTRACT 

It has been recently shown that molecular clouds do not exhibit a unique shape for 
the column density probability distribution function (A^-PDF). Instead, clouds without 
star formation seem to possess a lognormal distribution, while clouds with active 
star formation develope a power-law tail at high column densities. The lognormal 
behavior of the A^-PDF has been interpreted in terms of turbulent motions dominating 
the dynamics of the clouds, while the power-law behavior occurs when the cloud is 
dominated by gravity. In the present contribution we use thermally bi-stable numerical 
simulations of cloud formation and evolution to show that, indeed, these two regimes 
can be understood in terms of the formation and evolution of molecular clouds: a very 
narrow lognormal regime appears when the cloud is being assembled. However, as the 
global gravitational contraction occurs, the initial density fluctuations are enhanced, 
resulting, first, in a wider lognormal iV-PDF, and later, in a power-law A^-PDF. We 
thus suggest that the observed A^-PDF of molecular clouds are a manifestation of 
their global gravitationally contracting state. We also show that, contrary to recent 
suggestions, the exact value of the power-law slope is not unique, as it depends on the 
projection in which the cloud is being observed. 

Key words: ISM: general - clouds - kinematics and dynamics - turbulence - stars: 
formation 



1 INTRODUCTION 

It has been known since the first observations of molecular 
gas in star- forming clouds that the CO lines exhibit super- 
sonic line widths jWilson et al.l Il970l ) . iGoldreich fc KwanI 
l|l974l ) suggested that such supersonic linewidths could be 
produced b y large-scale collapse of th e molecular clouds. 
In contrast, IZuckerman fc EvansI (|l974l ) dismissed the idea 
of large-scale collapse with the argument that, if clouds 
were collapsing freely, the star formation rate in the Galaxy 
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should be a factor of 100 times larger than observed, lead- 
ing to quick exhaustion of its gas content. Those authors 
proposed instead that the large linewidths are produced by 
small-scale supersonic turbuIenc^U- 

Since then, the scenario of supersonic molecular 
cloud turbulence has received much attention, and tur- 
bulent models of molecular clou ds have b een devel- 
oped by several groups (see iMac Low fc Klessen. ,20o3 : 

^ Strictly speaking, IZuckerman fc EvansI l ll974h ruled out only 
radial large-scale collapse, which is hardly surprising, as molecular 
clouds are generally far from round. 
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IBallesteros-Paredes et al] |2007| . and references therein). 
However, the turbulent picture faces several problems, not 
the least of which is that supersonic shocks dissipate en- 
ergy quickly, within one dynamical time scale, even if tur- 
bulence is magnetized (IStone. Ostriker. fc Gammid 1 19981 : 
iMac Low erailil998l : IPadoan fc Nordlundlll999l l. Moreover, 
although stellar energy has been proposed as a source of 
turbulence within molecular clouds, it is not clear that such 
feedback can maintain the identity of the clouds; high-mass 
stars can easily blow apart clouds, and bipolar flows from 
low-mass stars are so focused that general support seems un- 
likely (see review bv I Vazq uez-Sem adcni 2010a). Thus it is far 
from clear - and it has not been demonstrated numerically 
- that stellar feedback can keep clouds near equilibrium for 
several dynamical timesca les. In fact, the available eviden ce 
suggest the opposite (see IVazguez-Semadeni et "all l20ld ). 

The idea that molecular clouds (MCs) are the result of 
large-scale compressions of the H I diffuse medium was ex- 
plored with numerical simulations of a 1 kpc squared piece of 
the galactic disk, representative of the Solar Neighb o rhood 
by [Ballesteros-Paredes, Vazquez- Semadeni & Scald (I1999I) 
and Balles toros-Parcd cs, Hartmann & Vazqucz-Scmadeni 
lll999h. and collected in a cohere nt sce nario by 
iHartmann. Ballesteros-Paredes fc" BerginI (120011 1. These 
authors found that MCs can be formed rapidly and pro- 
ceed to star formation almost at the same time, since 
the column density threshold for gravitational instability 
and f or molecular gas form a tion a re similai0. Simultane- 
ously, iHennebelle fc Peraultl (| 19991 ) discussed the idea of 
large-scale compressions of the diffuse medium undergoing 
thermal instability, thus producing dense gas. This led to 
the idea that molecular clouds in the solar neighborhood 
are formed in large scale compressio ns from the warm , 
diffuse, thermally unstable medium ([Hei tsch et al. ' '20061; 
IVazguez-Semadeni et akllioOTl : iHeitsch fc H artmann 2008). 
As these clouds accumulate mass and cool rapidly, they 
become Jeans unstable; at this point, gravity begins 
to dominate the motions, developing near equiparti- 
tion between the grav ita tional a nd kinetic energies 
JVazguez-Semadeni et akl l2007l : iHeitsch fc HartmannI 
I2OO8I ), errone ously but frequently infer red to be in "Virial 
Equilibrium" (|Ballesteros-Paredesll2006l ). As a consequence, 
the initial turbulent fluctuations in density and velocity 
produce hierarchical fragmentation, in which dense clumps 
with short free-fall times proceed to a rapid collapse, while 
at the same time the whole cloud contracts at a slower rate, 
because of its lower average density. 

This model of molecular cloud formation has been 
subject to d ifferent observational tests. For instance, 
iHeitsch et all (2009) showed that the CO(l-O) observed line 
profiles through such modeled clouds reproduce the super- 
sonic, turbulent line profiles of actual molecular clouds, 



^ Note, however, that the entire evolutionary sequence 
that culminates with the formation of a giant molecu- 
lar cloud involves timescales of a few tens of Myr, of 
which the " molecular stage" may represent a small fi- 
nal fraction jHartmann, Ballesteros-Paredes fc Ber ginI I2OOII : 
iBergin et ai]|2004l : IVazquez-Semadeni et al.ll2007l . l201of ) 



as well as the so-called core-to-core velocity dispersion. In 
other words, such global collapse of irregular structures 
develops internal disordered, turbulent motions at some 
level, and thus the large linewidths in actual MCs must 
be showing mainly the large-scale systematic, inhomoge- 
neous inward motions, ra t her th an pure turbulence. In fact, 
IVazguez-Semadeni et al] (|2009t ) showed that the culmina- 
tion of the collapse of a medium-sized cloud produces a 
dense clump with large velocity dispersion, with physical 
properties similar to those of massive star-forming regions. 
In this clump, the velocity field is dominated by the collapse 

motions. 

In addition, IBallesteros-Paredes et ahl ll201lf), based on 



observations of massive dense cores Jcaselli fc Mverd 



Plume et al] 1 19971: IShirlev et al] l2003l: ICibson et all 



1995; 



200E 



Wu et al.ir201Gl '). showed that the iL arsonl (119811 ) relation 



ship between velocity dispersion and size arises naturally if 
turbulent motions are the result of hierarchical and chaotic 
global -I- local collapse. In particular, they found that the 
velocity dispersion should scale as 



Sv 



2GT.R 



(1) 



where Sv is the velocity dispersion of the gas, E is the column 
density of the core, R its size, and G the uni versal constant of 
gravity, according to the Hove r et al.l (j2QQ9l 'l relation recently 
found for molecular clouds. 

In the present work we take a further step to compare 
models with observations. We analyze the column density 
probability distribution function from models of molecular 
cloud formation and evolution via warm, thermally unstable 
H I stream collision. The goal of the present contribution is 
to assess the evolution of the column density probability dis- 
tribution function (A'^-PDF) in these models, and compare 
it to the A''-PDF of observed molecular clouds. This is par- 
ticularly important because a number of recent observations 
have reported A^-PDF for several nearby molecular clouds, 
permitting more tests of theory. 

In section iJHwe discuss the current knowledge of volu- 
metric and column density pdfs for interstellar gas. In i]3.1l 
we briefly describe the models used, and in ^ we present 
the time evolution of our A'^-PDFs. In §3) we present a brief 
discussion and in iJSjwe give our main conclusions. 



2 VOLUMETRIC AND COLUMN DENSITY 
PDFS OF MOLECULAR CLOUDS 

One important tool for understanding the internal structure 
of molecular clouds is the volume density probability distri- 
bution (p-PDF). For isothermal, turbulent flows of a given 
Mach number M, where shocks occur stochastically and in 
succession, creating density fluctuations pi/po fx M^, the 
p-PDF is a lognormal functi on, i.e., a Gaussian functi on in 
the logarithm of the densitv. I Vazguez- Semadenil (|l994 ) and 
iPassot fc Vazguez-Semadeni] ~( 19981 ) explain this result as 
follows: in an isothermal flow, the speed of sound is spatially 
uniform, and a shock of intensity M will induce a density 
jump p\/ po over the mean density po. If another shock of in- 
tensity M arrives at the place where the density is now pi , it 
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will induce another density jump of size P2/P1 oc M^. Thus, 
if the flow is stochastic, such density jumps must be spa- 
tially uniform, and a given density distribution must be ob- 
tained by a succession of multiplicative density jumps, which 
are additive in the logarithm. Thus, at a given position in 
space, each density jump is independent of the previous and 
the following ones, and therefore, the central limit theorem, 
according to which the distribution of the sum of identically- 
distributed, independent events approaches a Gaussian, can 
be applied to the logarithm of the density, resulting in a 
lognormal density probability distribution function. 

Various authors suggested that the column density of 
turbulent, isothermal flows must als o be lognormal (e.g., 
lOstriker et al.ll200ll : [Ridge et al.ll2006l) . However, this is only 
a consequence of the fact that interstellar turbulence con- 
tains the largest velocities at the largest scales, and thus, 
the size of the interstellar clo uds is comparable to their cor- 
relatio n length. In this case, IVazauez-Semadeni &: Garcia I 
(|200ll ) showed that the local values of the volume density 
along a single line of sight are correlated, and thus the col- 
umn density might be representative of the mean value of the 
volume density along such line of sight. Thus, for different 
lines of sight, the values of the A'^-PDF are correlated to the 
mean values of the p-PDF, and thus, the A^-PDF follows the 
shape of the p-PDF. This would not be the case, however, 
if the correlation lengths of molecular clouds were small. 
In this case, the line of sight will contain many correlation 
lengths, and thus, the mean value will not be representative 
of the volumetric density of the line of sight. In this case, 
the central limit applies over the sum of several correlation 
length s, and the column density should ex hibit a Gaussian 
shape IVazguez-Semadeni fc Garcia Il200l} ). 

In recent years, a number of observational studies 
have focused on understanding the column density PDF of 
nearby clouds. For instance, Ridge_et al. (2006) used extinc- 
tion data from the COMPLETE (Coordinated Molecular 
Probe Line Extinction Thermal Emission) Survey of Star- 
Forming Regions, to argue that the gas in Ophiuchus and 
Pers eus molecular clouds exhibits a lognormal A'^-PDF(see 
also lGoodman et al.ll2009l ). Unfortunately, the column den- 
sity in this work is plotted in linear scale in the y axis 
(number of events), and thus, the excess in the high column 
density part of the A'^-PDF cannot be properly appreciated. 
However, that excess clearly indicates that some degree of 
departure of the lognormality is going on in Perseus. 

On the other hand, iFroebrich et al.l (|2007l ) presented 
maps and 7V-PDFs of 14 different regions inferred from ex- 
tinction measurements. They found that the A^-PDFs ex- 
hibit different shapes, although the dynamical range in vi- 
sual extinction Ay is frequently smaller than one order 
of magnitude, making it difficult to draw firm conclusions 
about the actual shape of the A''-PDF for these clouds. 

iKainulainen et ahl (|2009l ) took this a step further, show- 
ing that the Coalsack and Lupus V molecular clouds, which 
do not have signs of star formation, exhibit lognormal A'^- 
PDFs, while clouds that have already formed stars, like Tau- 
rus and Lupus I, do exhibit power-law like tails at large 
column densities . Their results were confirme d almost im- 
mediately afterwards by IFroebrich fc Rowlej (|2010l ). The 



standard interpretation of both groups is that clouds that 
are turbulent and non-star forming, exhibit a lognormal A'^- 
PDF, while clouds that are forming stars are somehow de- 
coupled, and exhibit a power-law tail at large column den- 
sities. 

The above observational results have a clear coun- 
terpart in numerical simulations of star formation 
in turbulent, isothermal MCs. Several groups have 
shown that the p-PDFs in those si mulation s de- 
velop power-law tails at late stages (iKlesseiil I2OO0I : 



iDib fc Bur kert"2005': 'Vazgu ez-Semadeni et all 1200^ . More 
recently, [Kritsuk et al. (201(j) have shown that centrally- 
peaked de nsity distri butions, such as the singular isother- 
mal sphere [Shul (|l977l ). or other dynamical solutions such as 
a pressure-free or inside-out collapse, which are all charac- 
terized by power-law density profiles, have p-PDFs and A*'- 
PDFs with power-law high-density tails, whose slope does 
in fact depend on the slope of the density profile. Also, they 
show that numerical simulations of self-gravitating isother- 
mal turbulence develop power-law tails as time progresses. 
Thus, they attribute the development of such tails in self- 
gravitating turbulent fiows to the formation of local collaps- 
ing sites. Although they do not mention it explicitly, one can 
conjecture that the PDFs of the entire turbulent flow exhibit 
the characteristic power-law tails of the collapsing centers, 
because the latter are the only sites where such high den- 
sities develop. Note also that, because those authors only 
evolved their simulations to less than half a free-fall time, 
the dominant component of the kinetic energy in the simu- 
lations was still that due to the turbulent flow. 

On the other hand, multiphase simulations of the 
formation and evolution of MCs out of the conver- 
gence of w arm, diffuse gas in th e presence of self- 

grav ity by IV azquez-Sema deni et al.l (|2007l . I2OIOI ) and 
Heit sch fc Hartm ann (2008) show that the dominant com- 
ponent of the kinetic energy in the clouds after they become 
dominated by self-gravity is generalized gravitational con- 
traction. In the remainder of this paper, we show that those 
multiphase simulations also develop power-law tails in their 
iV-PDFs at late times, suggesting that the PDFs of MCs 
should not be expected to be stationary, but rather to evolve 
in time, as a consequence of the transition of the clouds from 
being turbulence dominated to being c ollapse dominated. 

As an alternative interpretation, iTassis et al.l (|2010l ) 
have recently shown that a lognormal column density PDF 
does not necessarily imply that the flow is turbulent. 
In particular, they have shown that simulations of indi- 
vidual and multiple collapses also exhibit lognormal A'^- 
PDFs at early stages in their evolution. However, it is 
well known that isothermal c ompressible turbule n ce de - 
velops such a lognormal PD F Vazguez-Semadenil l|l994 ): 
iPassot fc Vazguez-Semadeni I (|l998l ). Thus, it is natural to 
ascribe the lognormal parts of PDFs of ISM turbulence sim- 
ulations to supersonic isothermal turbulence. Similarly, the 
high-density gas in the ISM is characterized by supersonic 
linewidths and behaves very close to isothermally, suggest- 
ing that lognormal PDFs in MCs are due to this process as 
well. 

In the next section we will show that this behavior can 
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be understood in a unified model of molecular cloud for- 
mation and evolution, where clouds are assembled by large- 
scale convergent flows in the WNM, which initially produce 
turbulence, but gravity rapidly takes over, producing a hi- 
erarchy of nested collapsing motions, and thus, the A^'-PDF 
develops a power-law wing at large column densities. 



3 RESULTS 

3.1 Brief description of the models 

In the present section we analyze the evolution of the N- 
PDF in numerical simulations of the formation and evo- 
lution of molecular clouds. In order to understand the 
evolution of the A''-PDF, we briefly describe the evolu- 
tion of the clouds of tho se simulations (for a review, see 
I Vazauez-SemadenillioiObI ) . 

The models consist in the collision of two large streams 
of warm thermally bi-stable H I gas. The streams collide 
with a given inflow velocity, on top of which fluctuations in 
magnitude and shape are superposed. 

As the streams collide, the compressed regions un- 
dergo a nonlinearly induced thermal instability, cooling 
down rapidly and becoming t urbulent due to a c ombi- 
nation of various instabilitie s (|Heitsch et al. I l2005l . l2006l : 
IVazguez-Semadeni et al]|2006l ). As the cold clouds accumu- 
late mass, they quickly become Jeans unstable at various 
scales. The superposition of those different centers of col- 
lapse produces hierarchy of nested col lapses, which is re- 
flected in highly supersonic line profiles l|Heitsch et ahllioogl : 
iBallesteros-Paredes et al]l201ll ) that traditionally have been 
interpreted as evidence of turbulence. 



VazQuez-Semadeni et al. 


20071 


Vazauez-Semadeni et al. 


l2010f) 



have performed simulations 



with different numerical schemes (SPH, fixed grid or AMR, 
respectively), and with slightly different initial conditions. 
However, in spite of the differences, the formation and evo- 
lution of the modeled molecular clouds exhibit similar char- 
acteristics: initially, all modeled clouds fragment rapidly due 
to dynamical an thermal instabilities, as well as the initial 
velocity fluctuations. However, while the models without 
gravity remain fragmented and do not develop large den- 
sity regions, the models with gravity start to collapse trans- 
versely to the inflow direction, developing local centers of 
collapse with large volumetric densities 10® cm~^), as 
well as large column densities, 1 0^"^ cm~^). 

As lHeitsch fc HartmannI (|2008l ) discuss, the same cloud 
at different evolutionary times appears as if it were differ- 
ent kind of objects: at earlier times, the compressed layer 
could be classified as a diffuse H I cloud, wit h low column 
densities and thus, very little CO con tent (see lHeitsch et al] 
l2006l : I Vazauez-Semadeni et all 120061 '). As time goes by, the 
layer accumulates more mass, column densities increase, and 
some CO begins to form. At this stage the cloud could be 
identified as a "translucent" cloud. Later in time, larger col- 
umn densities are achieved not only because of the compres- 
sion of gas along the direction of the streams, but because 



the transverse collapse causes a rapid increase of the col- 
umn density, shielding the interior of the cloud against UV 
photons. At this stage, the cloud appears mostly molecular, 
with dense cores that are already collapsing and which a few 
Myr later form stars. The precise times in which each one of 
these stages occurs depends on several parameters, e.g., the 
initial values of the density and velocity fields, the actual 
values of their fluctuations, the total mass of the streams, 
etc. However, once the cloud has accumulated enough mass 
to allow for molecule formation, the first collapsed objects 
appear a few Myr later. 

In the present work, in order to focus on the shape 
of the A'^-PDF during the process of cloud formation, we 
make use of runs without stellar fe edback. We use runs 
HF an d GF2 of the suite p resented bvlHeitsch fc HartmannI 
(|2008h and LAFO run by IVazguez-Semadeni et al.l l|2010h . 
Run HF is similar to GF2, but without self-gravity, and is 
used only for reference. Run GF2 represents the evolution 
of a ~ 22 X 44 X 44 pc cube filled with gas at T ~ 1800 K 
and at a density of n = 3 cm~^. The incoming streams are 
ellipsoidal cylinders with a cross section of 22 pc, an inflow 
velocity of 7.9 km/sec, and inflow boundary conditions on 
the numerical box, i.e., infinite incoming cylinders. In run 
LAFO, the simulation is performed using an adaptive mesh 
refine ment numerical scheme (ART code bv lKravtsov et akl 
Il997t ). The box has a size of 256 pc per side, filled with 
a gas at T = 5000 K, and a density of n = 1 cm~^. The 
streams are circular cylinders having a diameter of 32 pc, 
and an incoming velocity of 5.9 km/sec, and a length of 
110 pc. The inflows are thus fully contained within the box. 
This means not only that they have a finite length and 
duration, but also that t he to tal mass of the box in the 
IVazguez-Semadeni et al] (|2010|) runs is constant, while in 
the runs by Heitsch fc HartmannI (|2008l ) it is not. This re- 
sult will have implications on the height of the A'^-PDF, as 
we will see in tj3.2l 

In order to trigger instabilities in the compressed layer, 
in the HF and GF2 runs the shock front has a sinusoidal 
shape, while in LAFO run, a fiuctuating velocity field is 
added to the infiow velocity with rms amplitude of ~ 30% of 
the inflow speed, and with ~ 1/2 the diameter of the inflows. 
For more details, see conditions of runs GF2 and LAFO for 
iHeitsch fc HartmannI (120081 1 and IVazguez-Semadeni et al.l 
( 2010l ). respectively. 



3.2 Column density PDFs 

Figs. [1] [1 and [3] show the iV-PDFs of models HF, GF2 
and LAFO, respectively. For each run, we present two ex- 
treme cases: A'^-PDFs calculated by projecting the density 
field perpendicular to the infiow and one projecting along it. 
All plots exhibit 16 panels, from 3 to 15 Myr in Figs. [T] and 
(2] and from 5.6 Myr to 26.6 Myr in Fig. [3] In both cases, 
the time is measured from the beginning of the collision. The 
first point to notice in all three figures is that in all the simu- 
lations, at every time step, and in both projections, there is 
a large number of data points at low column densities. This 
is just a consequence of the setup of the simulations: the 
incoming fiows are cylindrical, and outside these cylinders 
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the gas is uniform, with low volume density and zero initial 
velocity. The peaks in each case simply reflect the column 
density of the uniform background medium. The mean den- 
sity of this field is ~ 1 cm"'', and thus, the typical column 
density of this part of the box is ~ 2 x 10^" cm"^ for run 
GF2 and ~ 8 X 10^° cm'^ for the LAFO one. 

The second point to notice is that the column densi- 
ties in both projections are extreme cases because these are 
extreme projections. Typically, molecular clouds must be 
observed at intermediate angles. In any case, the initial 7V- 
PDFs for the projection along the flow (plane y — z) axe nar- 
rower than the projections perpendicular to the flow (planes 
X — y ox X — z, not showed here). This is because the clouds 
are formed in the y — z plane, and hence the cloud is thinnest 
along the x-direction. 

In all three figures we further show lognormal and/or 
power-law fits as follows: Fig. [1] shows only lognormal 
fits (dotted line) at every timestep, since the correspond- 
ing run is purely hydrodynamical. Note that the PDF in- 
creases in height as time goes on. This is a consequence 
of the inflow bound ary conditions of the models from 
iHeitsch fc HartmannI (2008.' ), which imply that the total 
mass in the box increases in time. 

In Figs.[2]and|3l on the other hand, we show the lognor- 
mal fits from t = 5.3 Myr to t = 9.1 Myr and from t = 9.8 
Myr to t = 16.8 Myr, respectively. The final times of these 
intervals correspond to the time at which the cloud start 
collapsing more vigorously. However, for reference, we keep 
the first fit until the end of the simulation (dot-dashed line). 

We note that in both panels of Fig. [T]the A''-PDFs ex- 
hibit a lognormal-like wing at large column densities (dotted 
line) since early on in the simulation's evolution, and main- 
tain this shape for the rest of the evolution. The lognormal 
wing at large column densities is due to the density fluctu- 
ations produced by the turbulence induced by a combina- 
tion of nonlinear thin shell, thermal a nd Kelvin-Helmholt z 
instabilities in the compressed layer (|Heitsch et al.ll2005h . 
However, it is worth noting that this lognormal part is quite 
narrow, with column densities spanning less than two orders 
of magnitude, from ~ 2x 10'^^ cm~^ to 10^^ cm~^, for the en- 
tire duration of the evolution. This is due to the fact that the 
turbulence generated by the instabilit ies is not very strong, 
in ag reement with previous results ama fc Inutsuka I 
l200a ). 

The case with self-gravity is quite different (see Figs. [2] 
and [3}. We divide the behavior in two cases: the lognor- 
mal and the power-law behaviors. We first note that, after 
some transients, the dense column density region develops a 
narrow lognormal shape, indicated by the dotted line (which 
we mantain through all frames for clarity). However, as time 
goes on, gravity modifies the lognormal shape, first by mak- 
ing it wider. In order to show this widening more clearly, 
in Fig|4]we plot the evolution of the standard deviation a 
of the lognormal fits as a function of time for the four runs 
with self-gravity. From this figure it is clear that the lognor- 
mals tend to increase in width as the col lapse proceeds. Thi s 
behavior has been recently reported bv lTassis et al.l |20l3). 
Finally, at late times, gravity runs away and column den- 
sity grows several orders of magnitude in a few Myr. As a 



consequence, the A'^-PDF develops a power-law tail at large 
column densities. 

The A^-PDFs shown in Figs. [2] and [3] are remi- 
niscent o f the variety of A i'-PDF s for different clouds 
shown by Kainulainen et al.l |20l3). In fact, as shown by 
iFroebrich et al.l (|2007l ). every cloud may have its own slope, 
depending not only on the physical conditions of the cloud, 
but also on the resolution at which the cloud itself has been 
observed. 

Furthermore, by comparing the direction of the projec- 
tions, we note that the slope depends on the orientation of 
the cloud: although in models GF2 the slope in both projec- 
tions is almost the same (see Fig [5} , models LAFO (Fig [21 
show that the slope is steeper when we observe the cloud 
along its longer dimension. This effect occurs because if the 
cloud is mostly aligned with the line of sight, the observed 
column densities will be larger than the column densities 
obtained if the cloud is mostly on the plane of the sky. This 
result contrasts with the suggestion bv lKritsuk et al.l ()2010[ ) 
that all the A''-PDFs should tend to evolve towards a single 
value. These authors used a model of a spherical cloud, and 
thus every projection is the same. However, real clouds are 
far from round, and thus one cannot expect an asymptotic 
single value of the slope of the power-law tail of the A^'-PDF. 

Also, it is noteworth y that the other run with s lightly 
different perturbations in lHeitsch fc HartmannI l|2008l ). (run 
GFl, not shown), which has a smaller perturbation in the 
shock front, has an A^'-PDF evolution which is basically 
indistinguishable from the one presented here for model 
GF2. Similarly, the run with s maller amplitude ffuctua- 
tions (SAFO, also not shown) in IVazauez-Semadeni et al] 
(|2010h . follows an A'^-PDF evolution indistinguishable from 
model LAFO. However, the precise time in which each 
model ijHeitsch fc Hartmann|[200^ : IVazauez-Semadeni et al.l 
20101) start collapsin g are different: t ~ 10 Myr for the 
Heitsch fc Hartmann ||2008|) ru ns, and t ~ 15 Myr for the 



Vazguez-Semadeni et al.l 1 20100 runs. The difference seems 



to be due t o the free-fall time in each ru n: the compressed re- 
gion in the lHeitsch fc HartmannI (|2008l) runs has 47% of the 
volum e of the compressed region in Vazguez-Semadeni et al.l 
(|2010l ), but its mass flux is 5 times larger. Thus, the den- 
sity in the compress e d reg ion is 2.23 times larger in the 
iHeitsch fc HartmannI l|200^ run, which in turn implies a 
free-fall timescale 2/3 times smaller. This result suggests 
that, to first order, the main parameters that define the 
time of collapse are the velocity and th e density of the incom- 
ing field. Indeed, a parameter study bv lRosas-Guevara et al.l 
(2010) has found widely different star formation efficiencies 
and times of the onset of star formation as parameters such 
as the inffow velocity and the amplitude of the background 
turbulent fluctuations are varied. 

In both runs, after the development of the lognormal 
Ai'-PDF, a clear power-law tail develops at high column 
densities. This shape transition marks also the beginning 
of the accelerated large-scale collapse ongoing simultane- 
ously with the small-scale local collapses, and occurs in few 
Myr (^4 — 5 Myr), regardless of the particular setup of 
the runs. It seems thus that once the cloud has achieved 
enough shielding to form molecules, the cloud proceeds to 



© 2006 RAS, MNRAS 000, [TH9] 



6 Ballesteros-Paredes et al. 




10000 
1000 
100 
1 

1 

10000 
1000 
100 
10 

1 

„ 10000 
c 

o 1000 
100 
10 

1 

10000 
1 000 
100 
1 





3.G0Myr 








4.60f>/>r 




5. 30 My r ' 








[■),q;)K'yr 








fi,.-i;;Myr ■ 








9.9Gt>.'yr : 




:C.?OMyr 




"•.40Myr ' 




12.20Myr 




13.00K''yl : 

\ 








'I.SOMyr ■ 



10" 10^' 10" 10" 10^' 10" 10" 10" 10" 10" 10" IC 
Column density [cm"'] 



20,„21,„22,„23 ,(,20,021,022,023 , o20 , (,2 1 , ^22 , o23 , o20 , o" 1 Q^^ 1 O" 

Column density [cm"'] 



10"10" 1 0"1 



Figure 1. a) Time evolu tion of the column density probability distribution function of a run without self-gravity (run HF in 
iHeitsch fc Hartmannll2008l) seen perpendicular to the flow (edge on), b) Same as a), but for the projection along the flow (face on). 
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Figure 2. a) Time evolu tion of the column density probability distribution function of a run including self-gravity (run GF2 in 
iHeitsch fc Hartmannll2008h seen perpendicular to the flow (edge on), b) Same as a), but for the projection along the flow (face on). 

collap se in few Myr, as pointed out by IHeitsch fc HartmannI few clouds with no signs of star fo rmation (Coalsack, Lu- 
l|2008t) . The rapidity of this transition explains why most pus V, see iKainulainen et al"]|2009h . 
of the molecular clouds in the Solar Neighborhood exhibit 

young star forming regions (few Myrs), and there are only Thus, we conclude that as global-l- local col- 

lapses proceed, the number of high column density 
points in the A'^-PDFs increases, causing a high A''- 
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Figure 3. a) Time evolution of the column density probability distribution function of model LAFO in lVazguez-Semadeni et al 
seen perpendicular to the flow (edge on), b) Same as a), but for the projection along the flow (face on). 



4 DISCUSSION 
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Figure 4. Time evolution of the standard deviation of the log- 
normal function fitted in the two projections of the two runs with 
self-gravity. As can be seen, in all cases a increases with time, 
implying that gravity makes wider the lognormal of a purely tur- 
bulent gas. 



Franco fc Cox 



PDF tail. As pointed out by . 

Hartmann. Ballesteros-Paredes fc BerginI (|2001l . 

Heitsch fc HartmannI (|2008l ). in the Solar Neighborhood the 



19861 ): 

and 



pressure due to self-gravity becomes importat at about the 
column density necessary to shield the compressed region 
against UV radiation (10^^ cm~^), allowing the formation of 
molecules. This situation corr esponds to the 3rd row in th e 
rightmost column of Fig. 1 in lHeitsch fc HartmannI (|2008l ). 
and corresponds to the moment in which the cloud has 
already started to collapse in the direction perpendicular to 
the shock. At this moment, the densest cores of the cloud 
reach column densities of 10^^ cm~^. 



Froebrich et al.l l^2007^■ iKainulainen et all l|2009l ) and 
Froebrich fc Rowled l|201ol) . argue that the different shapes 



of the observed A'^-PDFs must be due to a change in the gov- 
erning physical processes i n molecular clouds. In particular, 
iFroebrich fc Rowlej l|2010t ) argue that most of their iV-PDFs 
can be fitted by a lognormal distribution, plus a power-law 
wing, the first one for low column densities, which corre- 
sponds to the gas that is not forming stars, and the second 
one for the gas which is forming stars. 

Our results show that these two regimes (quiescent 
vs. active star formation) may very wel l be part of a sin- 
gle ev olutionary process, as depicted by IKainulainen et ahl 
|200^. In the early stages of the molecular cloud, the 
H I colliding streams form the molecular gas and provide 
part of the supersonic turbulent kinetic energy observed 
in the line profiles. During this stage, the A'^-PDF is log- 
normal, and most of this stage occurs in the atomic phas e 
IVazquez-Semadeni et alll2007l : iHeitsch fc Hartmannll2008l ) , 
during which the cloud does not form stars. Later on, hi- 
erarchical and chaotic gravitational contraction takes over, 
tending to produce a power-law tail in the A^'-PDFs. The ini- 
tial manifestation of this effect is that the lognormal appears 
to become wider. It is important to mention that gravity is 
required in order to achieve the 1.5-2 orders of magnitude 
in the spread of the column density values, since pure tur- 
bulence produces only very narrow lognormal A''-PDFs. The 
contribution of self-gravity in increasing the column density 
helps the cloud to achieve values hi gh enough to rapidly 
form molecular gas, as poin ted out by lHeitsch fc HartmannI 
l|2008h : lHeitsch et all (|2009l ). A few Myr after the molecules 
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have formed, the densest cores collapse and the first stars 
begin to form. 

Although we a gree with the general scheme depicted by 
iKainulainen et all (2009), we stress that the contribution of 
self-gravity is indispensable to achieve high values of the 
column density. That is, self-gravity is not only important 
in driving the collapse of gravitationally unstable clumps, 
but also in forming such clumps, which woul d not be pro - 
duced by turbulence alone ( Vazquez- S em adcni et alll2008h . 
During the later stages, the primary source of "turbulent" 
motions (i.e., supersonic linewidths) in the molecular cloud 
is gravity as well, with local mass conce ntrations that render 
the velocity field more locally c omplex (|Heitsch et al.|[2"009l : 
iBallesteros-Paredes et"al]|201ll '). 



5 CONCLUSIONS 

Using numerical models of the formation and evolution of 
molecular cloud from the diffuse interstellar medium, we 
have confirmed the suggestion bv IKainulainen et al.l (|2009| ) 
that the A''-PDF transits from a lognormal shape at early 
times, when the kinetic energy in the clouds is dominated 
by their initial turbulence (produced by various instabilities 
in the compressed layer that becomes the cloud), to having 
a power-law tail at high A*' at late times, when their kinetic 
energy is dominated by gravitational contraction, which cul- 
minates with the formation of stars. Our results thus explain 
naturally recent observations showing that non-star-forming 
clouds exhibit lognormal A^-PDFs, while star-forming clouds 
exhibit power-law tails at high densities, understanding the 
result as the consequence of the transition of the clouds from 
a more diffuse, turbulence-dominated regime to a denser, 
star-forming, collapsing one. We also have shown that the 
slope of the power-law tail produced by gravity in the A'^- 
PDF does not have a unique value, but depends on the ori- 
entation of the cloud with respect to the line of sight. 
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